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SUMÁRIO 


Pretende-se no presente trabalho definir as bases de cálculo, isto é, a idealização 
da realidade que tem interesse prático considerar no projecto das barragens abóbada. 
Para isso, procura-se julgar a relevância da influência dos diversos parâmetros de que 
depende o estado de tensão deste tipo de estruturas. E discutida a influência das formas 
que podem ter as barragens abóbada, das propriedades do betão da estrutura, das 
características dos maciços rochosos de fundação e das solicitações actuantes, É tam- 
bém discutida a hipótese, correntemente aceite, das tensões normais se distribuirem 
linearmente ao longo da espessura das barragens abóbada. 


1. — INTRODUÇÃO 


Tratando este simpósio da teoria das barragens abóbada, compreende-se que os seus organi- 
zadores tivessem desejado que em primeiro lugar fosse formulado o problema físico envolvido, dada 
a sua grande complexidade. 

Vamos pois procurar nesta lição, por um lado, pôr em evidência os múltiplos parâmetros que 
influenciam o comportamento das barragens abóbada e, pór outro lado, definir a idealização da rea- 
lidade, isto é, o conjunto de hipóteses a partir do qual devem ser desenvolvidas as teorias de 
cálculo dessas estruturas. Tal como em toda a teoria física, essa idealização deve traduzir o que é 
relevante do ponto de vista da finalidade a atingir com a teoria, finalidade que, no caso presente, 
é o julgamento da segurança da barragem, o qual implica a determinação do estado de tensão desen- 
volvido. 

Os parâmetros que influenciam o estado de tensão das barragens abóbada, assim como de 
qualquer outra estrutura, são: 


— a forma da estrutura 
— as propriedades do material ou materiais constituintes 
— as ligações exteriores, isto é, no caso em consideração, as fundações 


* Versão em português da lição inaugural do «Simpósio sobre a teoria das barragens abóbada», organizado pela 
Universidade de Southampton, Inglaterra, em Abril de 1964. 
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— as solicitações actuantes, quer dizer, as acções de qualquer natureza que determinem defor- 
mações ou tensões. 


Além destes parâmetros os valores das tensões calculadas dependem ainda das hipóteses sim- 
plificadoras introduzidas no estabelecimento das próprias teorias de cálculo. 

Nos artigos seguintes vamos considerar sucessivamente todos os factores que acabam de ser 
mencionados. 


2. — AS FORMAS DAS BARRAGENS ABÓBADA 


As formas das barragens abóbada são muito variáveis, dependendo, por um lado, da forma do 
vale e das propriedades mecânicas das suas vertentes, e por outro lado, das concepções dos projec- 
tistas sobre o funcionamento estrutural e dos conhecimentos e meios de que eles dispõem para 
a elaboração dos projectos. Não há formas padrão, verificando-se mesmo uma tendência para a diver- 
sificação, em consequência da preocupação de busca de soluções mais económicas e do progresso dos 
conhecimentos e dos meios analíticos e experimentais de previsão do comportamento, donde resulta 
a possibilidade de uma adequação cada vez mais perfeita entre a estrutura e as condições do local 
onde vai ser construída. 

As concepções gerais sobre as formas e o funcionamento das barragens abóbada sofreram uma 
profunda evolução nas últimas décadas. Essa evolução ainda prossegue mas, segundo se nos afigura, 
de maneira mais atenuada, o que não é de estranhar se se atender a que as concepções cada vez 
mais se apoiam em conhecimentos de carácter científico e menos na intuição individual. Este factor 
está contribuindo para uma certa unidade de concepção geral que se sobrepõe à diversificação imposta 
pelas condições locais a considerar em cada caso. 

Nas concepções que primeiro surgiram, as barragens abóbada eram assimiladas a uma série 
de arcos sobrepostos, cuja interacção era desprezada, e ao paramento de montante era dada a forma 
cilíndrica, vertical, de directriz circular. Os arcos eram calculados pela fórmula dos tubos, donde 
resultavam barragens com arcos de espessura constante. Reconhecida a necessidade de ter em consi- 
deração o encastramento nas fundações, daí resultou a adopção de arcos de espessura variável, 
crescendo para as nascenças, em cujo cálculo era possível ter em consideração a deformabilidade das 
fundações. 

Um passo fundamental foi dado ao reconhecer-se o interesse, para a redução do volume da 

à obra, de dar aos arcos raios decrescendo com 

2 a cota, de harmonia com a diminuição da 
sen? 5 do Ha NNE | LI! abertura do vale. Daqui resultou uma modi- 
| DO 47 ficação profunda da forma, pois que se foi 

o “Ni obrigado a abandonar a forma cilíndrica ver- 


q 

Vs ad x Ed | tical do paramento de montante, aparecendo 

is 3 | assim a barragem de dupla curvatura. Cal- 
culados os arcos pela fórmula dos tubos, 

a conclui-se que, para um dado valor da ten- 
EE são desenvolvida, 7, e da corda do arco, c, o 

sm nEsam — volume mínimo de cada arco, V, e portanto 
” pd Rr nana da barragem, se obteria dando aos arcos um 
180 160 140 ma 100 80 8º ângulo ao centro de cerca de 133” (fig. 1), 


conclusão que muito tem influenciado as for- 
mas das barragens. Porém, a preocupação, 
que por vezes ainda se verifica, de adoptar 
tal valor não se justifica, De facto, a curva de variação do volume tem um mínimo muito pouco 
acentuado (fig. 1), de tal modo que, por exemplo, para um angulo 0 ==100º o acréscimo de volume 
é somente de cerca de 7 “, variação que não tem significado, dadas as grosseiras hipóteses em que 


Fig. 1 
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se apoia o cálculo do volume V. A tendência geral verificada nos últimos anos tem sido a de redu- 
zir os ângulos ao centro, com vista não só a obter soluções com menor volume de betão, quando 
dimensionadas por critérios mais evoluídos, mas também a melhorar as condições de inserção da 
barragem, isto é, a aproximar da normal à vertente o impulso descarregado sobre a fundação. Se se 
atender a que a segurança das barragens abóbada depende essencialmente das suas fundações, esta 
preocupação reveste a maior importância, devendo mesmo aceitar-se muitas vezes a sua satisfação 
à custa de incremento substancial das tensões desenvolvidas na estrutura. 

Na evolução das concepções sobre as formas e o funcionamento das barragens abóbada, outro 
marco fundamental é a consideração da real solidariedade que existe entre os arcos de uma barragem. 
Para dimensionar as barragens abóbada de harmonia com essa concepção desenvolveram-se os dois 
bem conhecidos tipos de métodos de cálculo, o designado por «trial-load» e os baseados na teoria 
das cascas, aos quais os computadores electrónicos trouxeram nos últimos anos apoio precioso. 
Simultâneamente foram desenvolvidas técnicas de estudo sobre modelo que são hoje correntemente 
usadas, sobretudo no projecto de obras importantes. 

A consideração da solidariedade entre os arcos tem permitido aproveitar melhor o material, 
isto é, obter soluções mais delgadas e mais económicas. Verifica-se um progressivo refinamento 
das formas, guiado pelo desejo de aproximação da forma ideal, que seria aquela para a qual a bar- 
ragem funcionasse como uma membrana submetida a uma compressão uniforme. O método expe- 
rimental da membrana fornece sugestões preciosas sobre as formas mais convenientes a adoptar 
para cada forma do vale (1). 

Para melhorar a distribuição das tensões constata-se a tendência nos últimos anos para adoptar 
arcos com raio de curvatura variável, crescendo para as nascenças, tendo sido usados arcos circu- 
lares de três, ou mesmo mais, centros, arcos parabólicos e arcos elípticos, em geral com espessura 
variável. Acresce que estas formas apresentam a vantagem de melhorar a inserção da barragem nas 
vertentes do vale. Deve notar-se que para vales muito estreitos, nos quais o efeito consola, isto 
é, o efeito da solidariedade entre os arcos, não é relevante, a adopção de arcos não circulares tem 
menos interesse. Também para vales com grande abertura parece ter 
pouco interesse a adopção de arcos não circulares, uma vez que o com- 
portamento dos arcos é sobretudo comandado pela deformabilidade das 
consolas (2). 

importante referir que a solidariedade entre os arcos acarreta, 
em certas zonas, o desenvolvimento de tracções sob a acção da pressão 
hidrostática, em especial na base a montante, na direcção vertical, 7pm 
(fig. 2), efeito este que é tanto mais acentuado quanto mais aberto é o 
vale. A preocupação de anular, ou pelo menos reduzir substancial- 
mente, estas tracções à custa do peso próprio constitui um factor deter- 
minante da tendência referida para uma certa unidade de concepção 
geral de formas. De facto, para conseguir tal resultado é actualmente 
muito corrente dar à consola central uma forma análoga à indicada na 
fig. 2, com um desaprumo para montante tal que o peso próprio ins- 
tala uma compressão 7ym suficientemente elevada para compensar a 
tracção devida à pressão hidrostática, sem que a jusante o peso instale 
uma tracção 74; importante. Além disso, é por vezes necessário dar um 
desaprumo para jusante na parte superior, a fim do peso próprio poder 
também compensar tensões de tracção verticais na zona A, instaladas 
pela pressão hidrostática. É oportuno referir que, apesar da concepção 
das barragens abóbada ser muito influenciada, sobretudo no caso de 
vales largos, pela preocupação de limitar o valor das tracções na base 
a montante, o facto é que a repercussão na segurança da estrutura da eventual fendilhação tia 
zona está longe de se encontrar esclarecida. 

Quanto às formas das barragens abóbada, é ainda digna de nota a tendência para aceitar vales 
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assimétricos, renunciando-se aos trabalhos de escavação necessários para assegurar a simetria. Tal 
é possível em virtude do desenvolvimento dos conhecimentos gerais sobre o funcionamento estru- 
tural das barragens abóbada e da capacidade dos meios de previsão do seu comportamento. No caso 
de vales com uma assimetria muito marcada é mesmo aconselhável que os próprios arcos da barra- 
gem sejam assimétricos e com curvatura mais acentuada na vizinhança da encosta mais inclinada, 
tal como é sugerido pelo método da membrana. 

Ao referir as formas das barragens abóbada a considerar nas teorias de previsão do seu com- 
portamento, não pode deixar de se mencionar as aberturas destinadas a descarregadores de super- 
fície, a descargas de meio fundo e de fundo, a condutas para alimentação de turbinas, etc. Nume- 
rosos ensaios sobre modelo de barragens abóbada com aberturas de variadas formas, dimensões e 
localização têm-nos mostrado que, dentro da precisão com interesse na prática, mesmo aberturas 
com dimensões importantes em face das da barragem, somente influenciam localmente o estado de 
tensão (3). Isto é, quando há aberturas em jogo elas podem ser ignoradas na determinação do estado 
geral de tensão da barragem, e o estado de tensão na vizinhança de uma abertura pode ser determi- 
nado isolando uma zona de barragem que a contenha e considerando as superfícies cortadas subme- 
tidas às tensões do estado geral. Porém, as concentrações de tensões em torno das aberturas consti- 
tuem em regra equilíbrios tridimensionais complexos cuja determinação só pode ser feita com a 
conveniente precisão, na maior parte dos casos, por via experimental. Quanto às tensões de compres- 
são desenvolvidas nestas concentrações, elas não levantam qualquer problema, o que é consequência 
dos valores relativamente baixos das tensões de compressão adoptadas no projecto das barragens 
abóbada, os quais permitem assim o desenvolvimento de elevados factores de concentração. Quanto 
às tensões de tracção, estudos diversos sobre modelo têm-nos evidenciado que as fissuras que pode- 
rão ocorrer não são em regra progressivas, limitando-se à zona onde se desenvolve a perturbação 
causada pela abertura, o que é compreensível, uma vez que o estado geral de tensão de uma barra- 
gem abóbada é de compressão. Portanto, as fendas que possam ocorrer não comprometerão a segu- 
rança da estrutura mas podem comprometer a estanquidade da obra. 

Outro aspecto a considerar é a frequente existência de encontros gravidade. Quando os encon- 


à 1 
tros têm uma altura importante em face da da barragem, digamos superior a 5 como ordem de 


grandeza, eles podem influenciar o estado geral de tensão da estrutura e portanto deve procurar-se 
dispor de métodos de cálculo com capacidade para os considerar. Mesmo para encontros com menor 
altura é necessário poder prever o estado de tensão na sua vizinhança nos próprios encontros. 

Um problema que por vezes se levanta é o da averiguação da influência, no estado de tensão 
da estrutura, de eventuais fendas que possam ocorrer em zonas onde se desenvolvam tensões de 
tracção, ou mesmo de juntas que se prevejam nessas zonas para nelas ter lugar a libertação das 
tracções. É por isso conveniente que os métodos de cálculo tenham possibilidade de tratar esse pro- 
blema. Tal como já se referiu em relação às aberturas, também diversos estudos sobre modelo nos 
têm mostrado que em geral a influência de fendas se faz sentir sômente na sua vizinhança (3). 

Uma questão da ordem geral que se põe ao definir o sólido barragem abóbada é a da eventual 
descontinuidade introduzida pelas juntas verticais entre blocos. De facto, em virtude da retracção do 
betão, devida essencialmente à libertação do calor de hidratação, e por vezes, também, de assenta- 
mentos diferenciais das fundações, essas juntas abrem, de tal modo que, apesar da cuidada injecção 
de que são objecto, se pode levantar a dúvida sobre a legitimidade de considerar uma barragem 
abóbada como um sólido contínuo, tanto mais que são raras as oportunidades de observar directa- 
mente a maneira como as caldas de injecção preenchem as juntas. A observação do comportamento 
das barragens abóbada, em especial as leituras de extensómetros colocados na vizinhança das juntas, 
revela porém que em regra as juntas não introduzem descontinuidade no campo de tensões. Um 
caso convincente que tivemos ocasião de estudar foi o de uma obra com cerca de 35 m de altura, 
com espessuras bastante reduzidas, cujas juntas, por se tratar de uma ensecadeira, foram injectadas 
com pouco cuidado, além de que também existiam juntas horizontais de betonagem fissuradas, con- 
forme evidenciam os repasses que se observam na fig. 3. Apesar da barragem mais se afigurar como 
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um conjunto de blocos ajustados, a medição de tensões e o seu confronto com resultados de estudos 
sobre modelos e de cálculos mostraram que a obra se comportava como um sólido contínuo (4). 


Fig. 3 


Refiramos finalmente o problema das formas das barragens abóbada no que respeita à deter- 
minação das tensões devidas ao peso próprio, solicitação esta que tem a particularidade de se ins- 
talar à medida que prossegue a construção da barragem, isto é, trata-se de uma solicitação que 
actua num sólido cuja forma varia continuamente. No caso dos blocos da barragem serem desa- 
prumados para montante, tal como na fig. 2, à medida que prossegue a construção eles tendem 
a deslocar-se para montante e a afastarem-se, podendo portanto o estado de tensão devido ao peso 
ser calculado por blocos verticais independentes. Quando a barragem é desaprumada para jusante, 
em toda a altura ou sômente na parte superior, há a tendência para os blocos se encostarem e assim 
constituirem uma abóbada que irá transferir lateralmente parte do efeito do peso, dependendo essa 
transferência do maior ou menor desaprumo da barragem, da abertura do vale e da maneira como 
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cresce a construção. Estudos recentes sobre modelo (5), mostraram que este efeito pode ser impor- 
tante, sobretudo no estado de tensão na parte inferior da barragem. 

Do que dissemos até agora sobre as formas das barragens abóbada, a conclusão de carácter 
geral a tirar é que os métodos de cálculo, para terem interesse para o projecto, têm de ser extrema. 
mente versáteis no que respeita à possibilidade de consideração de formas. De facto, uma barragem 
abóbada é em geral um sólido com raios de curvatura e espessuras variáveis de ponto para ponto 
e ainda com singularidades diversas, tais como aberturas, encontros gravidade, fendas e juntas. 
Além disso, na determinação do efeito do peso é por vezes necessário considerar a sua instalação 
progressiva numa abóbada de forma variável. 

oportuno frizar a carência que se verifica de estudos sistemáticos sobre a influência dos 
numerosos parâmetros característicos das formas das barragens abóbada. 


3, — PROPRIEDADES DO BETÃO 


Trataremos em primeiro lugar o problema das propriedades do betão a considerar nos métodos 
de cálculo destinados à previsão do estado de tensão sob a acção das solicitações previstas no pro- 
jecto. Vamos pois analisar em que medida é legítimo considerar o betão como um material homogéneo 
e obedecendo à lei de Hooke. Depois afloraremos o problema das propriedades do betão a conside- 
rar nos métodos de cálculo em relação à rotura. 

É bem sabido que, mesmo para valores muito baixos das tensões, o betão exibe deformações 
função do tempo. No entanto pode considerar-se, com suficiente aproximação, que até aos valores 
das tensões que se desenvolvem habitualmente nas barragens abóbada o betão é um material visco- 
-elástico e que as deformações devidas a sucessivas forças aplicadas são sobreponíveis (6). 

Sendo assim, pode demonstrar-se que se verifica uma consequência com grande importância 
prática: o estado de tensão de uma barragem de betão devido à aplicação de forças quaisquer pode 
ser calculado admitindo que o material obedece à lei de Hooke (7). Isto desde que se admita que 
existe um coeficiente de Poisson para as deformações visco-elásticas, o que pode ser admitido com 
suficiente aproximação. É oportuno referir a inconveniência de considerar nulo o coeficiente de 
Poisson, como é feito em certos métodos de cálculo, pois o seu valor tem uma influência apreciável 
no estado de tensão, sobretudo no valor das segundas tensões principais (8). 

Observa-se que, em rigor, o cálculo admitindo a elasticidade perfeita exige que as fundações 
exibam deformações função do tempo e que a relação entre os módulos de elasticidade do betão 
e da rocha se mantenha constante e igual à relação dos módulos de elasticidade instantâneas. 
Deve contudo notar-se que tal exigência não tem relevância, dado que, mesmo em rochas com um 
módulo de elasticidade bastante baixo — sômente as que influenciam apreciâavelmente o estado 
de tensão — a influência não é em regra importante, como se verá no artigo seguinte, sendo diminuta 
a sua repercussão no estado de tensão. 

Outra observação a fazer diz respeito à heterogeneidade da barragem resultante da diminuição 
da idade do betão da parte inferior para a superior e portanto da consequente variação das suas 
propriedades visco-elásticas. Daqui resulta um incremento da deformabilidade do material da bar. 
ragem de baixo para cima, que equivale de certo modo a uma redução progressiva de espessura 
Esta sômente poderá influenciar o dimensionamento da estrutura quando o betão da parte superior" 
for ainda muito jovem na altura em que é atingida pela água do reservatório uma cota vizinha 
da máxima, o que não é corrente, pelo menos em obras importantes. 

Um problema que se põe é o das propriedades do betão a considerar na previsão do efeito 
das variações de temperatura, as quais correspondem a uma imposição de deformações e não de 
forças. Então, a fluência do betão influencia em geral as tensões desenvolvidas, reduzindo-as, e 
tanto mais quanto mais lentamente se processar a variação de temperatura. Discutiremos adiante 
(art. 5) a legitimidade de substituir, como é prática corrente, a distribuição real de temperaturas ao 
longo da espessura das barragens abóbada por uma distribuição linear, desprezando o que é por 
vezes designado por efeito pele. Quer o efeito da distribuição linear de temperatura quer o efeito 
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pele serão influenciados pela fluência. No entanto, como as tensões desenvolvidas pelo primeiro são 
em regra pouco importantes em face das devidas à pressão hidrostática, o problema desse efeito 
não assume relevância. Porém, as importantes tensões instaladas pelo efeito pele podem ser substan- 
cialmente reduzidas pela fluência, para a onda anual de temperatura, como foi evidenciado por cál- 
culos de comportamento visco-elástico (9), mesmo no caso de betão com bastante idade. 

Em relação aos efeitos da temperatura, outro problema que se levanta é o da influência do 
calor de hidratação. Dada a elevada taxa de crescimento do módulo de elasticidade instantâneo e de 
decrescimento das deformações de fluência nas primeiras idades do betão, o efeito da fluência no 
estado de tensão provocado pelo calor de hidratação ainda é mais importante. Cálculos (10) e obser- 
vações têm mostrado que se dá mesmo a inversão do sinal das tensões que se desenvolveriam em 
regime elástico, podendo instalar-se tensões de valor bastante elevado. 

No cálculo do efeito do peso próprio o betão pode considerar-se elástico, uma vez que se trata 
da acção de uma força exterior. No entanto, no caso de se desenvolver durante a construção um 
efeito abóbada apreciável, pode acontecer que seja relevante a influência da heterogeneidade resul- 
tante da diferença de idades do betão, o que se traduzirá por uma transferência de tensões para a 
zona inferior da estrutura. 

Outra questão que se pode levantar é a da influência da heterogeneidade do betão no estado 
de tensão, em virtude da variação do seu módulo de elasticidade de ponto para ponto. Mesmo que 
se verifique um desvio quadrático médio de 20 */9 na resistência do betão — valor que não é prová- 
vel ser excedido — a ele corresponderá um desvio inferior a 10"/9 no módulo de elasticidade, o que 
significa ser improvável a ocorrência de desvios superiores a 20 "/o ou 30 “/o em relação ao módulo 
de elasticidade médio. Estudos em modelos heterogéneos mostraram que, mesmo variações do módulo 
de elasticidade como as que se acabam de referir, têm uma reduzida repercussão no estado geral de 
tensão. Se atentarmos a que os desvios considerados se referem aos ensaios correntes, sobre peças 
com volumes de poucos litros, e que, portanto, os desvios respeitantes a volumes significativos da 
estrutura serão muito mais pequenos, podemos concluir que não haverá necessidade de considerar a 
dispersão das propriedades no caso de toda a barragem ser construida com um mesmo tipo de betão. 

Finalmente consideremos o problema das propriedades do betão a ter em conta no cálculo das 
barragens abóbada em relação à rotura. Limitar-nos-emos a chamar a atenção para a influência da 

qualidade do betão nas deformações de rotura, 

y , R as quais têm um papel relevante no processo de 
RREO sacia e std redistribuição de tensões fora da fase elástica, como 
500 4 (é =1950210 tivemos ocasião de verificar experimentalmente. Na 


a o beto am j fig. 4 apresentam-se curvas obtidas em ensaios de 
E compressão de betões com diversas resistências (11). 

300 O = 210 kg/cm? Como se vê, a extensão correspondente à tensão 
” jE=2080x10' máxima é aproximadamente a mesma para todos cs 

= O = 89kg/cm 2 betões, mas a extensão de rotura aumenta substar- 
100 E - 1880 x 10º cialmente quando decresce a resistência do betão. 


Em resumo, quanto às propriedades do betão 
a considerar no projecto das barragens abóbada, 
pode afirmar-se ser legítima, nas condições habituais, 
a hipótese do betão ser um material homogénio e 
obedecendo à lei de Hooke no que respeita ao 
cálculo do efeito das forças exteriores. Quanto aos efeitos da temperatura, há necessidade, em 
certos casos, de considerar o comportamento viso-elástico do betão. 


0 1000 2000 3000 E (108) 
Fig. 4 


4. — CONDIÇÕES DE LIGAÇÃO EXTERIOR: FUNDAÇÕES 


A previsão do comportamento das barragens abóbada, como de qualquer outra estrutura, 
exige o conhecimento das condições de ligação exterior do sólido, isto é, no caso em consideração, 
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das relações entre as forças e os deslocamentos das fundações. Dada a acentuada heterogeneidade 
que em regra exibem os maciços rochosos, há em geral necessidade de considerar a deformabilidade 
como variável ao longo da fundação. Além disso, em cada ponto, as deformações são, na maior 
parte dos casos, função do tempo, da direcção de actuação das forças — consequência de aniso- 
tropia—e são muito influenciadas pela ocorrência de diaclases, falhas e outros acidentes geoló- 
gicos. Deste modo, a caracterização da deformabilidade das fundações constitui, muitas vezes, uma 
das maiores dificuldades que se apresentam ao projectista de uma barragem abóbada. 

Procuremos analisar qual o grau de influência da deformabilidade das fundações no compor- 
tamento das barragens abóbada, a fim de colher uma ideia geral sobre a «finesse» com que é neces- 
sário tratar o problema. Certamente que não é possível estabelecer regras rígidas, uma vez que essa 
influência depende de numerosos parâmetros, entre os quais citamos a extensão das zonas em que 
se verificam dados valores da deformabilidade, a forma do vale, especialmente a sua abertura 
e simetria, e a rigidez da barragem. Apesar destas limitações, procurámos num trabalho recente 
(12) sintetizar os resultados de numerosos estudos sobre modelo empreendidos pelo Laboratório 
Nacional de Engenharia Civil, de Lisboa, sobre a influência da deformabilidade das fundações. Tomando 
como situação de referência o estado de tensão de uma barragem abóbada quando a relação entre 
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os módulos de elasticidade da fundação e do betão, Es é tomada igual à unidade, concluímos, 
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necessariamente em termos muito gerais, que a variação do estado de tensão com a relação dos 
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mentalmente das fundações. 
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Fig. 5 


Uma primeira conclusão de ordem geral que se pode tirar da fig. 5 é que as barragens abó- 
bada não são tão sensíveis à deformabilidade das fundações como por vezes se pensa, podendo ser 
construidas sobre fundações com módulos de elasticidade muito baixos. Os resultados de estudos 
sobre modelos apresentados na fig. 6 confirmam a pequena influência da deformabilidade das fun- 
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as condições de fundação correntemente exigidas, variações de E, de 20"/y não têm em regra 
influência apreciável, isto é, bastará distinguir no cálculo da estrutura as zonas cujo módulo de elas- 
ticidade médio difere mais de 20 */o. 

Acerca da heterogeneidade das fundações não pode deixar de se frizar que em regra se podem 
observar deformabilidades muito diferentes em pontos próximos. É o que mostra a fig. 7 na qual, 
à esquerda, cada ponto tem por coordenadas os valores dos módulos de elasticidade determinados 
nas faces opostas, verticais, de galerias de ensaio, e à direita os valores no tecto e no pavimento. 
Os resultados referem-se a ensaios sobre uma superfície de carga de 1 m”, sendo de cerca de 4 m a 
distância média entre os bolbos de pressão interessados no par de ensaios correspondentes a cada 
ponto. Verifica-se que, para esta pequena distância, os módulos de elasticidade chegam a estar 
numa relação superior a 3. Porém, estas variações da deformabilidade têm uma influência muito 
localizada no estado de tensão da barragem e, dada a irregularidade da sua distribuição, nunca 
poderão ser consideradas no projecto. No caso de tais variações serem muito acentuadas poderá ser 
necessário dispor de armaduras ao longo da fundação, para absorção das tracções resultantes do 
efeito ponte sobre as zonas mais deformáveis, e assim evitar fendas, localizadas junto das funda- 
ções. Em conclusão, na prática só é possível considerar valores médios da deformabilidade de zonas 
de extensão apreciável da fundação. 

Não deve deixar de se frizar que na prática, mesmo com programas de ensaio «in situ» bastante 
desenvolvidos, é elevado o erro com que se pode garantir o valor do módulo de elasticidade de uma 
dada zona, excedendo na maioria dos casos 30'/o. Deste modo, é-se obrigado em geral a considerar 
valores de segurança, por defeito, para os módulos de elasticidade. O precário conhecimento dos 
valores reais dos módulos de elasticidade das fundações pode constituir, no caso de deformabilidades 
elevadas, uma causa de erro relevante na previsão do estado de tensão das barragens abóbada. 

Quanto à anisotropia das fundações, atrás referida, compreende-se, após o panorama que acaba 
de ser apresentado, que basta na prática considerar, em cada zona, a deformabilidade na direcção 
da força descarregada pela estrutura. 

Em relação ao efeito tempo na deformação dos maciços rochosos, a nossa experiência mostra- 

-nos que em regra não é relevante (12). Na 

EE (e) fig. 8 apresentam-se resultados de ensaios 
o «in situ», conduzidos nas condições já atrás 
referidas, designando-se por E e Eco respec- 
tivamente os módulos de elasticidade corres- 
pondentes a uma duração de aplicação de 
carga de 30 minutos e infinita. Como se vê, 
a diminuição relativa do valor de E é pe- 
quena, raras vezes ultrapassando 20/o. Con- 
clui-se, em face do que foi dito sobre o grau 


O 0,4 08 12 16 20 de influência da deformabilidade no estado 
E. 16º (kgem"?) de tensão das barragens abóbada, que em ge- 
ral não é necessário ter em consideração o 
ANTES DA INJECÇÃO DEPOIS DA INJECÇÃO efeito tempo. 
Quanto a falhas que atravessem a 
Granito o Granito o by Ê 
| | fundação, elas devem ser simplesmente enca- 
Xisto a Xisto : I 
radas como uma heterogeneidade do maciço 
Quartzite o Quartzite + 
: rochoso. Deste modo, e pondo de parte even- 
ig. 8 


tuais movimentos tectónicos, mesmo falhas 
com vários metros de caixa somente terão 
uma influência local. O mesmo se diz em relação às diaclases do maciço que pela sua importância 
devam ser consideradas à parte, isto é, cujo efeito não esteja integrado na'deformabilidade atribuída 
às diferentes zonas da fundação. 

Considerámos até agora o problema da influência da deformabilidade das fundações pondo de 
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parte a eventualidade da sua rotura. Porém, um problema de importância fundamental é o da averi- 
guação de segurança das fundações em relação à rotura tanto mais que, dentro das actuais concepções 
de projecto, a rotura de uma barragem abóbada só é concebível que possa ocorrer como consequência 
da rotura das fundações. Teria o maior interesse que se desenvolvessem métodos de cálculo em rela- 
ção à rotura que tivessem em consideração o funcionamento de conjunto da barragem e sua funda- 
ção, pois que na realidade elas constituem uma unidade estrutural. Presentemente, a previsão da 
segurança das fundações tem de ser feita a partir das forças descarregadas em regime elástico pela 
barragem, isto é, desprezando a redistribuição de tensões que ocorrerá no processo que conduz 
à rotura. Estudos sobre modelo mostraram que esta redistribuição é muito importante (13), sendo 
oportuno salientar que o método experimental tem possibilidade de ter em consideração as diaclases 
e falhas do maciço rochoso (14), as quais têm em regra uma importância decisiva na rotura. 

Acerca das fundações das barragens abóbada é oportuno frizar o contraste surpreendente que 
se verifica entre o cuidado exaustivo posto na previsão do comportamento e na execução da barragem 
própriamente dita e a pouca atenção dada à determinação das propriedades dos maciços de fundação, 
à consideração dessas propriedades no comportamento da estrutura e ao tratamento das fundações 
com vista a garantir as propriedades previstas no projecto (12). 


Para terminar este artigo queríamos destacar os seguintes pontos: 

— a deformabilidade das fundações tem em regra uma influência com interesse prático no 
estado de tensão para valores relativamente baixos do módulo de elasticidade das fundações; 

— quando estes valores ocorrem, a precária precisão do seu conhecimento pode comprometer 
a precisão de previsão do estado de tensão; 

— a heterogeneidade local que correntemente exibem os maciços rochosos, a qual não pode 
ser considerada nos cálculos, tem sómente uma influência local no estado de tensão; 

— interesse em que sejam desenvolvidos métodos de previsão do comportamento até à rotura 
da unidade estrutural barragem-fundação. 


5. — SOLICITAÇÕES : CONCEITO DE SEGURANÇA EM RELAÇÃO À ROTURA 


Quanto à definição das solicitações a ter em conta no cálculo das barragens abóbada, vamos 
considerar em primeiro lugar o caso do cálculo sob a acção das solicitações de projecto e depois o 
caso do cálculo em relação à rotura, observando-se que, actualmente, só é possível fazer o dimen- 
sionamento em relação a rotura por via experimental. 

O peso próprio de uma barragem e a pressão hidrostática, uma vez definido o nível da albu- 
feira, são solicitações sempre conhecidas com grande precisão. Quando se prevê o assoreamento da 
albufeira, já a definição da respectiva solicitação encontra a dificuldade de previsão do nível que 
será atingido pelo assoreamento e mesmo da distribuição das pressões correspondentes a um certo 
nível. Daqui podem resultar, em certos casos, indeterminações importantes para o estado de tensão 
da estrutura. 

Quanto à temperatura, são em regra conhecidas com a necessária aproximação as tempera- 
turas do ar, e também existem hoje elementos de informação satisfatórios sobre as temperaturas da 
água das albufeiras. Deste modo, e considerando a irradiação solar, pode ser calculada a distribuição 
de temperaturas, variável no tempo, ao longo da espessura das barragens abóbada. 

Já atrás foi referido que, segundo o critério correntemente adoptado, se substitui a distribuição 
real das temperaturas, t, ao longo da espessura da barragem pela chamada distribuição linear egui- 
valente, t' (fig. 9). Esta é tal que as deformações da estrutura são as mesmas que na distribuição 
real, uma vez admitida a hipótese da distribuição linear das extensões ao longo da espessura da 
barragem. No projecto das barragens abóbada desprezam-se pois as tensões devidas às variações de 
temperatura t-t', as quais não determinam senão um estado local de tensão, por vezes designado 
por efeito pele em virtude do campo de tensões assumir em geral sômente relevância junto da super- 
fície. Ora acontece que estas tensões atingem correntemente valores elevados (8), por vezes mesmo 
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ESPESSURA DA BARRAGEM 


Fig. 9 


superiores aos das tensões devidas à pressão 
hidrostática, quer de tracção quer até de com- 
pressão. Deste modo, não será de estranhar 
que não seja legítimo desprezar o efeito pele, 
mesmo contando com a plasticidade do betão. 
Ensaios realizados no nosso Laboratório so- 
bre peças prismáticas de betão, com gradian- 
tes térmicos ao longo da espessura, revelaram 
uma queda substancial de resistência à com- 
pressão (15) (Fig. 10). É de esperar que essa 
quebra de resistência se acentue quando cresce 
a resistência do betão, em virtude da conco- 
mitante redução da plasticidade. 

O problema ainda não está esclarecido, 
sendo possível que as tensões de pele influen- 
ciem a segurança das barragens abóbada em 
relação à rotura. A nosso ver trata-se de um 
dos problemas mais importantes a investigar 
sobre o comportamento das barragens abó- 
bada, sobretudo se se deseja aumentar o va- 


lor das tensões de segurança, como parece impor-se. À incerteza existente sobre o significado das 
tensões de origem térmica constitui uma das mais fortes razões para não ter interesse procurar-se 
determinar com elevada precisão as tensões devidas a outras solicitações. 

Uma questão que por vezes se põe é a necessidade da consideração da retracção higrométrica 
do betão no cálculo do campo de tensões das barragens. Dada a extrema lentidão com que se pro- 
cessa a secagem dos betões (16), para as espessuras correntes das barragens abóbada poderá, 
quando muito, verificar-se a secagem de uma camada superficial, dando origem a um efeito pele, 
Para as condições climáticas mais comuns este efeito poucas vezes interessará uma dezena de cen- 
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timetros. Acresce que a molhagem dos paramentos, a qual se propaga rápidamente ao interior do 
betão, pode eliminar o efeito pele em prazo muito mais curto do que o necessário à sua instalação. 
Deste modo, salvo muitos casos especiais, não há a considerar a retracção higrométrica do betão. 

Uma solicitação cuja consideração assume por vezes importância é a acção dos sismos, a qual 
tem sido assimilada à de forças estáticas definidas por critérios empíricos. Está hoje, porém, bem 
estabelecido o carácter de vibração aleatória que têm os sismos, impondo-se que sejam desenvol- 
vidos métodos de cálculo que o tenham em consideração. No entanto, uma das dificuldades com que 
ainda se luta é a própria definição, em termos de mecânica aleatória, da sismicidade de uma dada 
região, dada a carência de registos dos sismos de alta intensidade que interessam ao cálculo das 
estruturas. 

Considerámos até agora as solicitações actuantes na barragem própriamente dita, mas põe-se 
também o problema das solicitações actuantes directamente nas fundações. Entre estas a percolação 
da água da albufeira através do maciço rochoso pode ter particular importância na segurança das 
fundações em relação à rotura, em virtude não só das subpressões instaladas mas também da modi- 
ficação das propriedades do maciço rochoso (17). Quanto às deformações das margens e do fundo 
do vale devidas à pressão exercida pela água do reservatório, parece poder ter influência apreciável 
no estado de tensão das barragens abóbada a rotação sofrida pelas fundações (18). 

Quanto à determinação da segurança de uma barragem abóbada, como de qualquer outra 
estrutura, em relação à rotura levanta-se uma questão básica: como se deve supor que evoluem 
as diversas solicitações actuantes até se verificar a rotura? Este problema, de natureza muito 
complexa, apenas pode ser aqui aflorado. 

Em regra, uma estrutura está submetida a solicitações de diversa natureza física e, portanto, 
surge o problema de escolher as solicitações que vão variar e, para cada uma delas, se a sua inten- 
sidade variará uniformemente ao longo da estrutura ou não. O caso de uma barragem abóbada 
é muito representativo (13). As solicitações principais a considerar são o peso próprio e a pressão 
hidrostática. Conhece-se o peso próprio com muita precisão, isto é, o desvio em relação ao valor 
previsto nos cálculos é muito pequeno e assim não é pertinente considerar um coeficiente de segu- 
rança relativo a uma variação do peso. No que diz respeito à pressão hidrostática, não se justifica 
também, em regra, prever um aumento para além dos valores do projecto. No caso de uma albu- 
feira em que as encostas não sejam estáveis pode ser necessário considerar a possibilidade de escor- 
regamentos para a albufeira os quais provoquem vagas que galguem a barragem, aumentando assim 
substancialmente a pressão hidrostática. Foi o 


que aconteceu recentemente na barragem de dic A 
Vajont. Portanto se houver que recear tal a 
acidente, será conveniente determinar um coe- NÍVEL MÁXIMO / 
ficiente de segurança à rotura por aumento DE PROJECTO / , . 


da pressão hidrostática, definindo-se o valor 
desse coeficiente, por exemplo, como a razão 
a entre a pressão hidrostática máxima 
actuando na barragem quando ocorrer a ro- 
tura, h-, e a pressão hidrostática máxima 
prevista, h (fig. 11). Por vezes a segurança à 
rotura é avaliada supondo que a pressão 
hidrostática aumenta mantendo-se constante 
o nível de água na albufeira, definindo-se 
então o coeficiente de segurança como a rela- 


o Fr ao o 
ção Pr entre a pressão que determina a ro- 


tura, Pr, e a pressão hidrostática prevista, p | | á 1 
(fig. 11). Deve notar-se, contudo, que este Fig. 11 
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aumento da pressão hidrostática não tem contrapartida na realidade porque o peso específico da 
água não pode aumentar. Mesmo quando se prevê o assoreamento da albufeira, é ainda assim im- 
provável um aumento muito substancial da pressão para além dos valores do projecto. 

Outra possibilidade de rotura que não tem sido considerada (13) é consequência da resistência 
insuficiente dos materiais utilizados, que são, no presente caso, o betão e a rocha da fundação. 
Quando a consideração de um aumento das solicitações não tiver sentido, como acontece geralmente 
no caso das barragens, essa possibilidade é a única a considerar ao formular o problema de segu- 
rança à rotura. Então é possível definir um coeficiente de segurança à rotura como a relação entre a 
resistência prevista para o material da estrutura, 7, e o valor, 7%, que deverá ter para se dar a 
rotura permanecendo constantes as solicitações aplicadas. No caso das fundações em que a resistên- 
cia é em regra atrítica, pode tomar-se para coeficiente de segurança a relação entre o coeficiente de 
atrito esperado para o maciço rochoso, tg%, e o valor do coeficiente para o qual se dará a rotura, tg Or . 

Em resumo, na averiguação da segurança em relação à rotura devem ser definidos coeficientes 
de segurança de acordo com as condições físicas reais de funcionamento da barragem. 


6. — HIPÓTESES SOBRE A DISTRIBUIÇÃO DE TENSÕES 


Uma questão fundamental que se põe é a legitimidade da hipótese, admitida correntemente, da 
distribuição linear das deformações ao longo da espessura 
das barragens abóbada, hipótese que conduz a simplificações 
substanciais dos métodos de cálculo. Ou, por outras pala- 
vras, em que medida é necessário considerar as equações 
gerais a 3 dimensões da teoria da elasticidade. 

É de esperar que na vizinhança das fundações, muito 
especialmente na parte inferior das barragens e em encon- 
tros gravidade que eventualmente existam, a distribuição 
das deformações não seja linear, o mesmo acontecendo, cer- 
tamente, nas zonas singulares da estrutura, tais como aber- 
turas para descarregadores de superfície e para outros 
órgãos de descarga. O esclarecimento do problema luta com 
a dificuldade de medição de deformações em pontos do 
interior de modelos. Alguns resultados de estudos fotoelás- 
ticos a 3 dimensões, pelo método da congelação das tensões, 
confirmam a suposição da distribuição das tensões se afas- 
tar substancialmente da linear (18) na vizinhança das fun- 
dações (fig. 12). Na figura estão indicados os valores das 
tensões mormais devidas à pressão hidrostática, actuantes 
em secções horizontais da consola central de uma barragem 
com uma relação entre o desenvolvimento do coroamento e 
a altura de cerca de dois. 

A não linearidade das deformações ao longo da 
espessura junto das fundações deve ser particularmente 
acentuada no que respeita ao efeito da temperatura, pois 
há então ainda a contar com o fluxo tridimensional de calor' 


- TENSÕES DE COMPRESSÃO Dado o facto dos valores mais elevados das tensões 
+ TENSÕES DE TRACÇÃO nas barragens abóbada, quer de compressão quer de trac- 
Fig 12 ção, se desenvolverem a maior parte das vezes junto das 


fundações, ressalta bem claro o interesse da investigação 
do problema da distribuição das deformações ao longo da 
espessura nessa zona da estrutura. No que respeita à acção da pressão hidrostática, como se verifica 
um acordo satisfatório, para fins práticos, entre os valores das tensões fornecidas pelos métodos de 
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cálculo mais evoluídos e das medidas em pontos da superfície de modelos, conclui-se que a hipótese 
da linearidade não é muito grosseira, pelo menos no respeitante ao campo de tensões na superfície. 
Exceptuamos desta afirmação a zona inferior das barragens, sobretudo com maiores espessuras, onde 
temos observado discrepâncias importantes, os encontros gravidade e as zonas singulares strás 
mencionadas. O comportamento da zona inferior é particularmente importante pois, como já foi 


referido, é frequente desenvolverem-se nela tracções que condicionam a própria concepção da 
estrutura. 


7. — CONCLUSÃO 


Apesar de nos termos limitado a referir sumáriamente os parâmetros que influenciam o com- 
portamento das barragens arco, julga-se ter ficado bem marcada a extrema complexidade que pode 
assumir o problema da previsão do estado de tensão destas estruturas quando se queira considerar 
todos esses parâmetros. Põe-se porém uma questão de base, a saber, a precisão com que tem sen- 
tido na prática fazer essa previsão, dadas as incertezas envolvidas e as elevadas margens de segu- 
rança exigidas. Apesar de se reconhecer a dificuldade desta questão, mesmo de a analisar, vai-se 
procurar dar-lhe uma resposta, tanto mais que só assim será possível fazer um julgamento de valor 
das teorias de cálculo à disposição dos projectistas. 

Os erros de que vêm eivados os valores das tensões — isto é, os desvios entre os valores 
realmente desenvolvidos na estrutura em serviço e os valores calculados — dependem, por um lado, 
do esquema da realidade adoptado, quere dizer, das hipóteses definidoras do sistema físico em consi- 
deração e, por outro lado, dos erros introduzidos pelos métodos usados para o cálculo das tensões, 
incluidos nestes erros os resultantes das hipóteses preliminares sobre a distribuição de deformações 
ou tensões. 

Do que foi dito nesta lição depreende-se que os erros resultantes da esquematização da reali- 
dade que tem de se adoptar na prática podem ser bastante elevados, sobretudo os provenientes 
da definição das fundações e das solicitações térmicas. Sendo assim não se justifica procurar métodos 
de cálculo de elevada precisão. O que é muito importante é dispor de métodos com grande versati- 
lidade, isto é, que permitam tratar a extrema diversidade de situações que ocorrem na prática, em 
especial no que respeita às formas das barragens abóbada e às condições de fundação. Quanto à pre- 
cisão a exigir aos métodos de cálculo, acresce ainda que, por razões complexas que não serão anali- 
sadas, se exigem margens de segurança muito elevadas e definidas de maneira pouco precisa, mesmo 
arbitrária, o que também favorece a aceitação de métodos de menor precisão. 

Deve no entanto notar-se que os valores das tensões calculadas têm sobretudo um significado 
relativo, para confronto de diversas situações que se encarem para um mesmo local e mesmo de obras 
diferentes. Deste modo, e uma vez que se admita que os erros introduzidos pela esquematização 
da realidade são da mesma ordem — o que se verificará em especial no caso do confronto de solu- 
ções para um mesmo local — já se justifica uma maior exigência para a precisão dos métodos 
de cálculo. Tudo considerado, julga-se que, no estado actual dos conhecimentos sobre o comporta- 
mento das barragens abóbada, é aconselhável exigir, no caso de obras de certa importância, que 
os métodos de cálculo introduzam erros que não excedam cerca de 10 e 20"/o, respectivamente, nos 
valores das maiores tensões de compressão e de tracção devidas à actuação de todas as solicitações. 
A apreciação da precisão de um método de cálculo exige estudos sobre modelos reproduzindo diversas 
situações típicas, sendo oportuno referir o grande interesse do programa de investigação lançado 
pela Civil Engineering Research Association com vista a apreciar os métodos de cálculo. Deve fri- 
zar-se que somente em situações simples conseguem os métodos de cálculo correntemente usados 
alcançar o nível de precisão que acabamos de sugerir. 

Para terminar desejamos somente acentuar o grande interesse que têm as barragens abóbada 
do ponto de vista da economia e da segurança dos aproveitamentos hidráulicos e a possibilidade da 
sua utilização em vales muitos largos, sob a forma de obras de abóbadas múltiplas de grandes 
vãos (19), o que justifica todos os esforços para serem desenvolvidas teorias de previsão do seu 
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comportamento que permitam considerar todos os parâmetros relevantes e forneçam os resultados 
com a precisão conveniente. 


17. 
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SUMMARY 


The aim of the present work is to define the bases of calculation, i. e. the idealiza- 
tion of reality which presents practical interest to consider in the design of arch dams, 
For this purpose it is attempted to judge the relevancy of the influence of the varied 
parameters on which the state of stress in thistype ol structures depends. The influence 
is discussed of the shapes the arch dams may have, of the properties of the concrete 
used in the structure, of the characteristics of the foundation rock masses and of 
the acting loads. lt is also discussed the hypothesis, currently accepted, of the normal 
stresses linear distribution along the thickness in arch dams. 
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TRAÇADO DE DIAGRAMAS DE RENDIMENTO 
EM COLINA PARA TURBINAS HIDRÁULICAS 


ror JOSÉ REBELLO DA SILVA 


Engenheiro de Máquinas (I. S. T.) 
Bolseiro da Fundação Calouste Gulbenkian 
Equiparado a Bolseiro pelo 1. de Alta Cultura 


RESUMO 


Este artigo estuda a obtenção dos diagramas de funcionamento das turbinas defi- 
nindo as grandezas que nele intervém. Refere-se em particular 20 caso das turbinas 


Kaplan, 


1— INTRODUÇÃO 


O diagrama em colina de uma turbina hidráu- 
lica traduz as suas condições de funcionamento 
em todos os pontos de serviço. 

Este diagrama compreende essencialmente o 
traçado das curvas de igual rendimento num sis- 
tema de eixos coordenados n':, Q', podendo, 
adicionalmente, ser completado com as curvas 
de igual potência e com as de igual binário; no 
caso das turbinas Pelton podem ainda dese- 
nhar-se as curvas de abertura constante do in- 
jector ; no caso das Francis, além das curvas de 
abertura constante do distribuidor, podem dese- 
nhar-se, para as de rotação específica elevada, as 
curvas de valor constante do coeficiente de cavi- 
tação de Thoma; na caso das Kaplan, além des- 
tas duas famílias de curvas ainda se podem 
desenhar as curvas de abertura constante das 
pás da roda. (!) 

Os diagramas são determinados em laborató- 
rio sobre modelos reduzidos e são referidos 
normalmente a uma turbina unitária, isto é, uma 
turbina semelhante ao protótipo mas com o 
diâmetro de referência (*) D==1 m e funcionando 


(1!) Normalmente na zona de exploração efecfiva de 
uma turbina Kaplan existe uma relação bem determinada 
entre as aberturas do distribuidor e da roda e não é então 
necessário mais do que o traçado de uma das correspon- 
dentes famílias de curvas (ver também a nota de pé 6). 

(?) O diâmetro de referência varia com a firma cons- 
trutora. Como mais apropriado, o que não significa que 


sob a altura de queda H==1m, Esta turbina é 
caracterizada pelas grandezas unitárias: 


— velocidade: ny = ti (1) 
VH 
— caudal: Qu = lc (2) 
D'.VH 
— potência: Pq = E er cas (3) 
D2.VHº 
M 
— binário: M'|= ——— 4 
in I  H (4) 


que se deduzem das grandezas correspondentes 
do protótipo, de diâmetro D e funcionando sob 
uma altura de queda H, introduzindo as rela- 
ções de semelhança. Nas grandezas unitárias 
o índice 1 indica a referência a 1 m de altura de 
queda da turbina e a plica a referência a 1 m 
de diâmetro. 

O traçado dos diagramas para as grandezas 
unitárias (embora no laboratório se não trabalhe 
na realidade, como regra geral, com turbinas 
modelo unitárias, mas sim com modelos de diá- 
metro bastante inferior a 1 m e funcionando sob 
alturas de queda variáveis e normalmente maio- 
res do que 1 m) tem a vantagem de permitir 


seja o mais utilizado, indica-se: para as turbinas Pelton, 
o diâmetro do jacto; para as Francis, o diâmetro de saída 
das pás da roda; e para as Kaplan, o diâmetro da roda. 
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uma passagem imediata ao diagrama do protó- 
tipo pela introdução das relações (1) a (4). Como, 
porém, o número de Reynolds do protótipo é 
superior ao do modelo correspondente, as perdas 
por atrito são proporcionalmente maiores neste 
do que naquele; ao traçar o diagrama da turbina 
protótipo é necessário, por isso, introduzir o 
respectivo «ganho em rendimento» (*) devido à 
sua maior dimensão. Observa-se, no entanto, 
que o traçado dos diagramas das turbinas unitá- 
rias é geralmente efectuado directamente com os 
valores dos rendimentos que foram determinados 
para a turbina modelo, sem a consideração do 
respectivo «ganho em rendimento». 

Do diagrama em colina do protótipo podem 
obter-se não só as curvas características Q,n e 
P,7 para valores constantes de H, mas também 
todas as restantes curvas e diagramas necessários 
ao estudo completo do funcionamento da uni- 
dade. 


2 — REALIZAÇÃO DOS ENSAIOS 


O equipamento de ensaio para a determina- 
ção de rendimentos compreende fundamental- 
mente as seguintes partes: instalação da tur- 
bina modelo própriamente dita (espiral, distri- 
buidor, roda, veio e tubo de aspiração, para as 
turbinas Francis e Kaplan; conduta, injector, 
roda, veio, carcaça e poço de descarga, para as 
turbinas Pelton); dispositivos para medição do 
binário (mecânico, hidráulico ou eléctrico), do cau- 
dal (diafragma, venturímetro, descarregador aferi- 
do, etc.), da altura de queda da turbina (manóme- 
tros, flutuadores, etc.) e da velocidade (taqui- 
metros electrónicos ou electro-mecânicos) e dis- 
positivos acessórios para completar o circuito 
(bombas, condutas, depósitos, canais, etc.). 

O que se segue refere-se a ensaios no caso 
de uma turbina Kaplan — já que a partir dele se 
pode deduzir a aplicação às turbinas Francis e 
Pelton —e compreende apenas o traçado das 
curvas de igual rendimento. 

O traçado destas curvas num diagrama em 


(') O «ganho em rendimento» — designação que atri- 
buímos à correcção dos valores dos rendimentos devido 
à maior dimensão da turbina protótipo em relação à tur- 
bina modelo - é determinado por meio de fórmulas apro- 
priadas, das quais as mais correntemente utilizadas são 
as de Ackeret, Moody e Hutton. 
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colina é efectuado a partir do conhecimento das 
curvas Kaplan Q%,'7%i;, correspondentes a dife- 
rentes valores de n'. Por sua vez, para o tra- 
çado de cada uma das curvas Kaplan é neces- 
sário o conhecimento das respectivas curvas 
hélice. Os passos a considerar são, portanto, 
pela ordem que se referem: traçado das curvas 
hélice, a partir dos dados obtidos nos ensaios; 
traçado das curvas Kaplan, a partir das curvas 
hélices; traçado do diagrama em colina, a partir 
das curvas Kaplan. 

A realização dos ensaios para a determinação 
das curvas hélices pode ser conduzida segundo 
dois métodos distintos. 

No primeiro método — que, de um modo geral, 
é o mais indicado para instalações de ensaio em 
circuito fechado — fixam-se, para cada turbina 
hélice, isto é, para cada posição estabelecida para 
o ângulo 7 de abertura das pás da roda, sucessivos 
valores convenientes da velocidade unitária n', 
e para cada um destes variam-se as posições de 
abertura ao das pás do distribuidor. Para cada 
conjunto de valores 7,n'1,ao, determinam-se, 
com base nas indicações dos dispositivos de me- 
dição as seguintes grandezas: 


— altura de queda da turbina, H (m), 

— caudal da turbina, Q (m'º/s) 

— velocidade da turbina, n (r.p.m.), 

— binário interno da turbina, Mi (m kg). (') 


Com estas grandezas calculam-se : 


— a poência interna da turbina, 


P = w), 
102 
— a potência da turbina, Bota (kW), 
102 
— o rendimento interno da turbina, 
102 . Fi 
nj = 
RS | 
dei o o 
— o caudal unitário, Q1= =— 
D2y H 


(*) O apoio do veio da turbina e do freio é feito de tal 
modo —em geral com a utilização de chumaceiras osci- 
lantes — que o valor medido no freio corresponde ao biná- 
rio interno da turbina: Mj = M + Mm em que M é o bi- 
nário da turbina e Mm o valor correspondente às perdas 
externas. 


No segundo método — que, em geral, é o mais 
conveniente para instalações de ensaio em cir- 
cuito aberto — fixam-se, para cada turbina hé- 
lice, sucessivas posições de abertura das pás do 
distribuidor e varia-se, para cada uma destas, o 
valor da carga no freio. Para cada conjunto de 
valores 7, ao, carga no freio, calculam-se, como 


5 E , 
anteriormente, os valores de P;, 4 Se pts 4 


calculando-se ainda o valor da velocidade uni- 
tária, n4 = ca : 
vH 

A diferença essencial entre as duas técnicas 
é que, enquanto na primeira se determinam para 
cada turbina hélice as condições de funciona- 
mento, para cada valor constante de n', fazendo 
variar os valores de ao, na segunda as condições 


ce funcionamento para cada turbina hélice são 


a; V, =Ps 


dentes a cada valor de n', a partir das grande- 
zas calculadas. 

O mesmo não acontece quando os ensaios 
são realizados de acordo com o segundo método, 
sendo então necessário um traçado preliminar 
das curvas n4, Qu en, 1; para se proceder 
em seguida ao traçado das curvas hélice Qi, 1. 

Considera-se o caso correspondente à realiza- 
ção de ensaios de acordo com o segundo método 
assinalado. 

Para cada turbina hélice, quer dizer, para cada 
valor 7x fixado para abertura das pás da roda, 
traçam-se as curvas n'1, Q'4 relativas a cada um 
dos valores (ao)y = constante, atribuídos à aber- 
tura das pás do distribuidor ; tem-se, pois, x dia- 
gramas e em cada um y curvas (*). Em sobrepo- 
sição traçam-se as curvas n',7; relativas tam- 
bém a cada um dos y valores de ao. 


(ae) 
(ao); 
(Quss|--------—=--4-=——-— (ao 
(A Ds PERNA SS > s4 (ao)2 
(Asa > (da), 
(a)sap----> E sis 
(o)sip o <L AR => 
“a [Mt%) 
RS a, Pe fi 
o = pa EN 
E RR 
o) h )3 
bo 
0 E meio: MA 
Fig. 1 


determinadas fazendo variar o valor de n', para 
cada valor constante de ao. 


3 — TRAÇADO DO DIAGRAMA DE RENDI- 
MENTO EM COLINA 


A realização dos ensaios segundo a técnica 
indicada como primeiro método permite o tra- 
çado directo das curvas hélice Q', ni; correspon- 


Na figura 1 representam-se as curvas n', Qi 
e n'4, Yi para cinco posições diferentes de aber- 
tura das pás do distribuidor e relativas a um 
determinado valor 7; de abertura das pás da 


eme . ——— 


(*) Esclarece-se que o número de curvas pode ser dife- 
rente para cada um dos x diagramas e que, além disso, 
aquelas curvas não correspondem necessáriamente a va- 
lores iguais atribuídos a a, . 
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roda. A partir destas curvas determinam-se os 
valores correspondentes de Q', e 1; para valores 
constantes de n'; na figura estão indicados os 
valores correspondentes de (Q')sy e (nilsy para 
valores constantes A e B de n4 (para não tornar 
o desenho confuso, os valores de (Q'1)y e (ni)sy 
correspondentes ao valor B de n', foram apenas 
assinalados sobre as curvas). 


(0), (Qpm 


(Qu 


origem dos diagramas traçados com os diferentes. 
valores de » sobre a mesma vertical, pois sendo. 
idêntica a escala em todos aqueles diagramas, a 
determinação dos valores de Qi e 7; para cada 
valor constante de n' é feita sempre sobre a 
mesma vertical de abcissa igual ao valor de 
nt. 

Na figura 2 traçaram-se as curvas hélice para 


(Quiv CQp)y 


( Qus: 4 Ass ( Q%)ss 
(O )s; (Q Os 


Fig. 2 


A partir do conjunto de todos estes diagramas 
traçam-se, para cada determinado valor cons- 
tante atribuído a n'4, as x curvas hélices Qu,r. 
Facilita-se esta passagem com o alinhamento da 


um valor A de n', e assinalaram-se sobre a curva 
“3 Os pontos A; a Às; em correspondência com 
os dados da figura 1. 

A envolvente de todas as curvas hélice repre- 


Fig 3 EN b 
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senta a curva Kaplan para o valor de n', consi- 
derado. 

O traçado do diagrama em colina é efectuado 
agora a partir das curvas Kaplan dos diagramas 
anteriores ("): 

Na figura 3 representam-se algumas curvas de 
rendimento constante, indicando-se os pontos | 
a V relativos ao valor À de n' e obtidos em 
correspondência com os dados da figura 2. 
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SUMMARY 


The author discusses the methods of obtaining the characteristic curves of a tur- 
bine, defining the interested variables. It is also taken in account the specific case of 


the Kaplan turbines. 


A OBSERVAÇÃO DA ONDULAÇÃO 


CARACTERÍSTICAS DOS TRENS DE ONDAS OBSERVADOS NA FIGUEIRA DA FOZ 


ERRATA 
Página Linha Onde se lê: Deve-se ler: 
464 2 da nota (!) resistência existência 
474 figura 14 m/sec? m/seg* 
476 18 de relação da relação 
477 17 Hiro H no 
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NOTAS INFORMATIVAS G. D. U. 624.344.5/94 
Elementos sobre a produção e o consumo de energia 
na rede eléctrica nacional 
— Elementos extraídos das estatísticas mensais do Repartidor Nacional de Cargas (R.N.C.) 


Nora: As produções e os consumos das empresas do R, N. C. representam 
cerca de 94 */, dos totais do Pais. 
I —Breve nota mensal 


Depois da ocorrência, em Março, de um regime de 
afluências superior à média — o que conduziu a um 
armazenamento no conjunto das albufeiras superior ao ABRIL 
respectivo valor da curva guia de segurança — regres- 
sou-se, em Abril, a um regime de afluências desfavo- 
rável o que levou a que se iniciasse, no dia 21 deste 


mês, o desarmanezamento do conjunto das albufeiras. 2.3.- DIAGRAMAS DE CARGA DOS DIAS CARACTERÍSTICOS PyI-EM-Ch, 
As afluências ao conjunto do sistema decresceram | 

rápidamente durante o mês de Abril por efeito da es- nas 

cassez de precipitações, notando-se, contudo, uma su- + 

bida transitória das afluências na passagem de Abril - + A. 


para Maio, originada pelas precipitações então ocorri- 
das. 


= 
; A 
Il — Elementos gerais (GWh) A f | 


Acumulados desde 1 de Janeiro 4 HA. 


+ Variação u SETE 
1964 | 1965 | op 
! 
Produção hidráulica (Ph) .. ./1529,1 13664 — 
Produção térmica (P+)... bo 4| 142,8 + 140 
Produção total (PT). ..... 1588,5 1509, 9i— 65 
Euergia recebida de empresas | y 
não pertencentes ao RNC (Er) | 14,9) 129 — 13 TE res | 
Exportações (Ex) .. ......l| 55,4 10, — 9 
Importações (1) , «« «casas 00 60.2 — 
Saldo importador (S1) ...... —36,4 +59,2 — 
Consumo em bomba gem (C JC!) — | 6,8 | — 
Produção para con- (4) 


sumos perman. (Pcop) (2) 1249, O 14169. + 13,4 
Produção para con- | : 
sumos não perman. (Pep) (?) 3180. 157,6 | — 0 
o, PA 15670 15819 + 1 
Coeficiente de hidrauticidade | LN 0,82 | E 


IV— Energia armazenada nas principais albufeiras 


No fim do mês 
Albufeiras: RE and 


NOTAS GWh 0% (1) 
(1) Energia consumida em bombagem para armazenamento na Al- o = n e 
bufeira do Alto Rabagão, | | 
(2) A produção para consumos permanentes (Pop) & determinada Alto Rabagão A AR | 215,5 221 
pela seguinte expressão ; Po = Pr4- E, + Sj — Cg== Penp Daradiio | 1104 496 
(4) Consideram-se consumos não permanentes os indicados em 4.2. E Ney de BE | ' , 
(4) O aumento percentual de produção para consumos permanen- Venda NOVA «vv co ko 115,0 sos 
a Sa ad atinge dias especiais, é Salamonde ar te a edpnadoa 22.8 | 892,6 
DO o a na 27,7 | 83,1 
III — Diagramas de carga dos dias característicos Ve sivmssics.s 48,11 46,5 
| 4: teira: Cabril A LR: DO ERR o CR 331,1 97,6 
| 1964 | 1965 Castelo do Bode. . .,..| 1252 76,8 
= Pa GUGA sms nasc cú TA 89,2 
hr pi 4% um 12661 | 14056 Lagoa Comprida . +... 29,5 (2) 67,5 
rodução térmica | b s LESS ss & os. o. ã 
Produção total (PT) MWh| 12651 | 14056 ro ja 90,1 
Trocas com [ Export. (Ex) MWh| o 0 ERNNRE id bio) io oie | Ef dt) 76,7 
Espanha Import. (1) ) MWh | O Õ POVOR sa via icms sa 6,5 (3) 49,0 
“Consum em bomb. hidroel. (Cp) MWh o TETRA Total com A, Rabagão . . .| 1092,6 51,8 
Prod, para cons. perm. (Peo) MWh| 9859 | 11388 sem A, Rabagão . . .| 877,1 77,2 
Prod, para cons. não perm, m. (Pero) MW 2762 | 2568 
TOTAL PT + (I-Ex) MwA | 12661 | 13956 - ani 
nua Potência máx, o | 683 | 7 . 
2 s| Pr + (Ex) Etdania min. MW, 3585 37 Notas: 
117 5 
: a = pe ed cm rs a | | Eco (!) Coeficiente de enchimento em energia das albufeiras, 
EE Potênciamax. MW| 562 “43 | (2) Inclui 3,3 GWh armazenados em Vale do Rossim no 
e á Potência min. 235 | 969 início do mês e 2,7 GWh no fm do mês, 
Um ni Utiliz. da ponta nr 17,5 17,7 | (º) Inclui 1,4 GWh armazenados no açude do Poio no 
Factor de carga | 0,73 0,74 


início do mês e 1,4 GWh no fim do mês, 


ANTIGAMENTE : 


HOJE : 


Com 
imagem 
direita 

na luneta 


Trabalho cómodo e simples 
sem erros e sem dúvidas, 
sinais correctos, não inver- 
tidos, claros e rápidos. 


O pequeno nível WILD N 10 com a 
designação «E» tem imagem direita do 
campo de observação na luneta. À mesma 
tem luminosidade igual ao modelo con- 
vencional e o mesmo tamanho. 


Representantes exclusivos 
WILD PORTUGAL, LDA. (tissoaz- rELEFONE cria 
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berliet-tramagal 


PARA TODAS AS CARGAS, PARA TODAS AS ESTRADAS 


MENALUILGICADUANRTE FERKRELEHIFLA,S. AR... 
divisao berliet RUA TOMÁS RIBEIRO, 50-A -LISBOA e RUA PASSOS MANUEL 36-PORTO 
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C. D. U. 628.392 


PRÉ-DIMENSIONAMENTO DUM EXUTOR SUBMARINO 
DE ESGOTO DE ÁGUAS RESIDUAIS 


ror AMÍLCAR JOSÉ RAMOS AMBRÓSIO 
Eng.º Civil (1.5. T.) ' 


do Gabinete de Estudos e Projectos 
de Engenharia Sanitária e Hidráulica 


O Engenheiro Armando Coutinho de Lencastre, de quem o autor deste artigo é 
colaborador no Gabinete acima referido, foi incumbido, pelo Governo Geral de Angola, 
da elaboração do Plano Geral de Saneamento de Luanda. 

Entre os vários problemas que houve que se resolver, um consistiu na fixação do 
esquema de saneamento do esgoto de águas residuais que, quer do ponto de vista 
técnico como económico, melhor se adaptasse aos condicionamentos existentes. 

Para tanto analisaram-se diversas soluções, uma das quais previa o lançamento no 
mar do esgoto proveniente duma extensão considerável da Cidade, através dum exutor 
submarino enraizado num determinado local da Ilha de Luanda. 

A fim de se computar os encargos correspondentes a tal solução, designadamente 
as parcelas relativas a investimento e custos de energia de bombagem, teve de se pro- 
ceder ao dimensionamento da obra o qual, como se verá, não podia deixar de ter um 
carácter aproximado. 

Uma vez determinado o diâmetro do exutor, estimado o respectivo comprimento 
e fixadas as demais característicss (terminal, tipo de material, revestimentos, protecção 
na zona de rebentação, etc.), foi possível fazer um cálculo suficientemente aproximado 
desses encargos. 

Verificou-se que, mesmo tendo em conta a margem de imprecisão resultante do 
dimensionamento aproximado feito, a solução não se apresentava econômicamente com- 
petitiva com a que finalmente veio a ser preconizada, 

Não obstante ter sido posta de parte, levou-se a sua análise, como aliás a de todas 
as outras soluções admitidas, a um grande pormenor uma vez que tal procedimento se 
considerou indispensável para se poder estabelecer conscienciosamente o esquema de 
saneamento mais aconselhável para a Cidade. 

O presente artigo tem por base o estudo então feito, limitando-se, todavia, aos 
aspectos relacionados com a determinação do comprimento do exutor e com a forma 


do terminal. 


1. INTRODUÇÃO 


O lançamento do esgoto no mar por exutores 
é, hoje em dia, pelo menos nos Estados Unidos 
da América, quer na costa atlântica como na costa 
do Pacífico, prática relativamente corrente, de 
tal modo que só na Califórnia, cerca de 125 ci- 
dades litorais, entre as quais se contam 11 das 
13 maiores cidades daquele Estado, recorrem a 
tal processo de rejeição e depuração natural do 
esgoto. [1] 

A técnica tem, também neste capítulo, evoluído 
bastante nos últimos anos. Alguns dados poderão 
sublinhar esta afirmação. 


Assim, dentro da maior variabilidade de com- 
primentos e de diâmetros dos exutores, de pro- 
fundidades finais a que são lançados e de descar- 
gas diárias que permitem, foram já atingidos os 
seguintes máximos: [1) 

a) — de comprimento: 11 500m num dos exu- 

tores de Los Angeles; 
b) — de diâmetro: 3,60 m num outro exutor de 
Los Angeles; 

c) — de profundidade: 100 m na Cidade de 
Everett, no Estado do Oregon; 

d) — de descarga máxima: 1 600 000 mº/dia de 
esgoto, pelo exutor de Los Angeles indi- 
cado em b). 
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a ——— LIMITE DA ZONA SERVIDA PELO EXUTOR 
ARA Edil 0 05 1 Km 


Fig. 1 — Esquema geral da solução analisada 
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De notar que o exutor de Los Angeles com 
3,60 m de diâmetro e dimensionado para 
1600000 m/dia, tem cerca de 8000 m de ex- 
tensão e vai até 60 m de profundidade. 

Deste modo, a situação a resolver em Luanda, 
que deveria prever o lançamento diário de 
28000 mº de esgoto por um exutor de 0,75 m 
de diâmetro com uma extensão oscilando em 
torno dos 1300 m, a uma profundidade de cerca 
de 40 m, não constituiria, do ponto de vista 
estritamente técnico, um impasse. A recorrência 
a uma empresa especializada seria o caminho 
para a resolução do problema. 


* * 


Na solução analisada o esgoto é bombado 
para o exutor tendo-se previsto, para o efeito, a 
construção de duas estações elevatórias e dum 
conjunto de condutas como, esquemáticamente, 
se apresenta na Fig. 1. 

Por conveniente localização das estações e tra- 
çado das redes de colectores, procurou-se tirar o 
máximo partido das condições topográficas exis- 
tentes e desse modo fazer afluir por gravidade 
as estações elevatórias referidas o esgoto da 
maior área possível da Cidade. 

Ao esgoto das restantes zonas fixou-se destino 
diferente daquele, não tendo interesse, todavia, 
no âmbito deste artigo, fazer-lhe qualquer refe- 
rência. 

O sistema de saneamento adoptado no Plano 
Geral é o do tipo separativo; nestas circunstân- 
cias, o esgoto que ocorre nas estações elevató- 
rias é apenas o de águas residuais. Dado o tipo 
de ocupação urbanística da zona interessada no 
exutor (predominantemente residencial e de ser- 
viços), o esgoto de águas residuais a lançar no 
mar era, na sua quase totalidade, esgoto domés- 
tico e como tal pode ser considerado. 

Antes de ser bombado para o exutor, o 
efluente seria sujeito a trituração por dispositi- 
vos apropriados a instalar nas estações elevatórias. 

A população servida pelo sistema descrito 
atingia 123 500 habitantes, correspondendo-lhe 
um caudal de ponta, debitado pelo conjunto dos 
grupos elevatórios das duas estações, de 730 1/s. 


* k x 


O enraizamento do exutor foi previsto na zona 
mais proeminente da Ilha, como se indica nas 
Hgs. 1 e 7. 


O lançamento do esgoto no mar do lado exte- 
rior à Baía pareceu evidente. Com efeito, a alter- 
nativa de lançamento na Baía não se mostrou 
viável na medida em que as profundidades sus- 
ceptíveis de provocarem uma boa difusão inicial 
do esgoto (superiores a 20 m) se encontram a 
partir de grande distância (cerca de 1700 m) da 
estação elevatória mais próxima. Para além disso, 
a intensidade e direcção de correntes, os fluxos 
e refluxos das marés e a forma fechada da Baia, 
dificilmente permitiriam que, mesmo com lança- 
mento a distância superior à mínima indicada, 
ficasse garantida a ausência de poluição nas 
zonas destinadas a banhos e a desportos náu- 
ticos. 

Assente o lançamento do esgoto do lado exte- 
rior à Baía, outra variante se poderia encarar, a 
qual consistiria em lançar o esgoto não por um 
único exutor mas por dois — um por cada estação 
elevatória. 

Ao primeiro, referente à EEI, corresponderia 
o local já indicado para a solução estudada. 

O segundo, dada a situação da EE2, teria o 
seu enraizamento frente a esta estação elevató- 
ria. Todavia, como se verificou, para que o ponto 
de lançamento do esgoto se situasse a uma dis- 
tância da linha da praia definida pelo prolonga- 
mento da Ilha para sul, sensivelmente igual à 
determinada na solução que se analisou em por- 
menor, seria necessário dar a tal exutor um com- 
primento superior a 2500 m a partir do seu 
enraizamento junto à EE2. 

Tal variante não se apresentou competitiva 
com a solução que foi objecto de estudo, dado 
os custos de construção dos exutores. 


2. DADOS OCEANOGRÁFICOS 
DISPONÍVEIS 


Não se dispunha de elementos suficientes para 
se fixar, com a máxima garantia, o comprimento 
conveniente para o exutor; eles eram, todavia, 
bastantes para se estabelecer uma gama de com- 
primentos dentro da qual a extensão real do 
exutor se fixaria. 

O processo corrente de dimensionamento de 
obras desta natureza compreende uma primeira 
fase de trabalho de gabinete, a única que então 
se podia realizar, destinada à estimativa do com- 
primento do exutor, à determinação do seu diá- 
metro e à definição das restantes características. 
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Se as conclusões a que se chegasse provassem 
que a solução era econômicamente viável, não 
obstante a margem de erro implícita no dimen- 
sionamento sumário, teria então que se levar a 
efeito uma série de trabalhos laboratoriais e de 
campo, que no caso mais geral deveriam englo- 
bar os seguintes aspectos : [2,3] 


a) — exigências requeridas, do ponto de vista 
de concentração bacteriológica, turvação, 
matérias flutuantes, depósitos, redução do 
teor de oxigénio, etc, às águas do mar 
destinadas a banhos; 

b) — análise do esgoto bruto, nomeadamente da 
sua concentração bacteriológica, B.O.D,, 
teor em matérias sólidas, etc; 

c) — ensaios in natura destinados a deduzir o 
poder bacteriofágico das águas do mar, 
isto é, da sua capacidade de redução da 
concentração bacteriológica e ainda des- 
tinados a medir o oxigénio dissolvido, 
tempo de regeneração, etc; 

d) — medição de salinidades e avaliação de 
temperaturas nas zonas interessadas na 
depuração do esgoto, quer à superfície 
do mar, como a várias profundidades, ao 
longo dum período de meses conveniente ; 

e) — ensaios com matérias corantes ou radioac- 
tivas, para avaliação da dispersão do es- 
goto e consequente diluição ; 

f) — análise de correntes, ondulação, forma- 
ção de vagas, etc; 

£) — levantamento batimétrico na extensão pre- 
vista para o exutor e respectiva sondagem 
dos fundos para conhecimento da sua 
natureza. 


O simples enunciado dos trabalhos a realizar 
deixa prever os meios de que se deveria dispor, 
quer em pessoal qualificado como em aparelha- 
gem, e concluir daí que efectivamente só empre- 
sas especializadas os poderiam levar a efeito. 

No tocante a dados oceanográficos foi muito 
útil a consulta que se fez do trabalho publicado 
pela Comissão de Estudos da Erosão e Defesa 
da Restinga de Luanda [4], comissão constituída 
em 1953 com o objectivo de determinar a viabi- 
lidade e a forma de suster a erosão da res- 
tinga pelo mar e traçar um programa de estudos 
capaz de levar ao completo conhecimento do 
fenômeno e à adopção de soluções definitivas. 
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Para cumprir a incumbência, a citada Comis- 
são teve que arrolar toda uma série de elementos 
existentes relativos a batimetria, marés, ventos, 
etc, proceder a uma série de observações de 
ondulação, marés, correntes, sedimentologia, etc, 
e fez executar um levantamento hidrográfico, 
traçados de perfis transversais e longitudinais, 
fotografias aéreas, etc, de toda a zona interes- 
sada no problema em causa. 

De todo o material existente e compilado e de 
todos os elementos posteriormente obtidos por 
iniciativa da Comissão, destacam-se, aqui, com 
interesse imediato para a questão do exutor, as 
observações de correntes. 

Para estudo do regime de correntes na zona 
vizinha da costa exterior da Ilha, realizaram-se 
seis tipos de observações, algumas das quais 
meramente qualitativas. 

De todos eles, o que compreendeu a observa- 
ção de correntes por meio de flutuadores em 
quatro verticais, é o que, pela sua própria natu- 
reza, pelo número de observações efectuadas e pelos re- 
sultados obtidos, é mais digno de menção [4]. 

Assinalam-se na Fig. 7 as quatro verticais onde 
se efectuaram as observações de correntes. Como 
se poderá verificar, são as verticais Vi, V2 e V3 
as que mais interesse têm para o estudo do 
problema do exutor, dado que se situam na 
zona onde, por um lado, ele poderia ser lan- 
çado e, por outro, onde se daria a dispersão da 
mancha de esgoto. 

A vertical Vi estava situada à profundidade 
de —4me a cerca de 250 m da linha da praia; 
aV:a —4m, também, e a cerca de 400 m; para 
a Vs aqueles valores eram de — 10 m e de, apro- 
ximadamente, 1000 m. 

O tempo de observações decorreu entre 21/4 
e 30/6 de 1954 e entre 6/9 do mesmo ano e 
3/1 de 1955. O período de trabalho em cada dia 
ia das 7 às 16 horas, lançando-se flutuadores 
com intervalos, sensivelmente regulares, de 
20 minutos. 

O número de observações realizadas foi de: 


894 na vertical Vi 
912 na vertical Vs 
1002 na vertical Vs 


ou seja, um total de 2808 nestas três verticais. 
Estes números mostram a amplitude do traba- 
lho levado a efeito, 


As observações em cada vertical foram efectua- 
das não apenas à superfície mas a várias pro- 
fundidades. 

Consultadas as cadernetas onde se fizeram os 
registos das observações e compilados os ele- 
mentos referentes a velocidades médias e má- 
ximas e número de observações, para vários 
azimutes, mas três verticais citadas, estabele- 
ceu-se a seguinte série de valores a utilizar em 
ulteriores cálculos: 


azimutes velocidades 
O sucos“ COMA 
O assuma zs 80 q 
OO” emos HD dá 
O” ss ums q x 06. & 
A” qn Rs os q ê 90 » 
160” avessas es DO =» 


OO sw as é 200 q» 
BOM da via. .+» 805 q 
MO ac rmus sa 8) & 


Admitiu-se, por hipótese, que, para cada azi- 
mute, a corrente era rectilínea e de intensidade 
constante e igual ao valor respectivo registado 
na relação anterior. 


3. LANÇAMENTO DO ESGOTO NO MAR 
POR EXUTOR 


3.1 — Índices de poluição 


Está provado, pela experiência que existe, de- 
signadamente nos Estados Unidos, dos problemas 
relacionados com o lançamento do esgoto no 
mar, que se se provocar uma boa difusão inicial 
do esgoto na massa líquida que circunda as saídas 
dos exutores, a qualidade da água nas zonas de 
banhos passa a ser exclusivamente governada 
pelo teor da concentração bacteriológica. [5] 

Será pois através do controle do teor da con- 
centração bacteriológica que será definida a ex- 
tensão do exutor. 


* x 


As bactérias que têm sido aceites, até agora, 
como melhor definindo o grau de poluição das 
águas, são as do grupo coliforme (B. coli), parte 
das quais, as chamadas bactérias fecais (Esche- 
richia coli), vivem no trato intestinal do homem. 
Deste modo, a determinação do grau de poluição 
duma água é, no fundo, circunscrita à contagem 


do teor de B. coli. Recentemente, porém, tem-se 
verificado não ser o grupo de bactérias colifor- 
mes o indicador bacteriológico de poluição mais 
seguro, nomeadamente nas águas do mar, mas 
somente a espécie Escherichia coli desse grupo; e 
há quem indique, em alternativa, o Streptococcus 
faecalis do grupo das bactérias enterococos, para 
além de outras várias bactérias que têm sido 
propostas como indicadores de poluição. As ra- 
zões destas divergências radicam-se, fundamen- 
talmente, no facto de, estatisticamente, se obser- 
var uma mais estreita correlação entre os teores 
de E. Coli e de S. faecalis e o número de doentes 
de carácter gastero-intestinal, entre os banhistas, 
ao mesmo tempo que se verifica uma grande 
variabilidade no teor de B. Coli. [6| 

Para além dos métodos laboratoriais de detec- 
ção de batérias em amostras de águas, a sua 
contagem faz-se, modernamete, por recorrência 
ao número mais provável (NMP) que representa a 
densidade de organismos que, em termos proba- 
bilísticos, existe muma determinada amostra. 
Greenwood e Yule, desenvolveram a teoria do 
NMP, e formularam as expressões matemáticas 
necessárias à sua determinação. [7] Determinados 
os NMP de várias amostras sucessivas, o teor da 
concentração bacteriológica da água em estudo 
é depois expresso em termos da média aritmética, 
ou da média geométrica, ou da mediana, etc, 
dos vários NMP calculados em cada amostra. 

Nos Estados Unidos, cada um dos Estados tem 
as suas próprias normas que estipulam os teores 
de concentração bactereológica em ordem aos 
quais se pode considerar aceitável a água do 
mar destinada a banhos. A variabilidade dessas 
normas levou a que se fizesse um inquérito des- 
tinado aos Departamentos de Saúde dos vários 
Estados, no sentido de, para além do mais, se 
averiguar em que bases se apoiavam as normas 
em vigor quanto à fixação dos teores de concen- 
tração bactereológica. [8] Invariâvelmente as res- 
postas a esta questão indicaram que na fixação 
desses teores não havia qualquer fundamento de 
ordem científica ou técnica, mas apenas o facto 
de a elas terem correspondido, até à data, con- 
dições de sanidade satisfatórias entre os banhistas. 
As normas dos vários Estados, referidas umas 
em termos de média aritmética, outras em média 
geométrica, etc, do NMP de B. coli; mostraram 
que, quando reduzidas ao mesmo índice (média 
geométrica do NMP de B. coli), as suas exigên- 
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cias se escalonavam entre 2,5 B. coliijml (Cali- 
fórnia) e 30 B. coli/ml (Michigan). 

Todas estas considerações, àcerca do índice de 
poluição bacteriológica, do modo de expressão 
dos teores de concentração e das normas que 
definem as características aceitáveis para as 
águas do mar destinadas a banhos, têm por fina- 
lidade mostrar a fluidez e mesmo as incógnitas 
que ainda existem em torno deste problema, 
apesar dos avanços da técnica, pelo que, no 
estudo feito, se considerou razoável supor que a 
média geométrica de 30 B. coli/ml nas águas do 
mar é índice de condições aceitáveis para banhos. 


3.2 — Concentração bacteriológica no esgoto 


O número de B. coli nas feses humanas é 
muito grande: a excreção diária, por pessoa, 
pode ir de 125><10º a 400><10º B. coli, rela- 
cionando-se os valores mais elevados com os 
meses mais quentes. 

Para a capitação média ponderada de 200 1/ 
/hab/dia, que se tomou no trabalho, se suposesse 
uma dispersão uniforme dos B. coli no esgoto, 
ter-se-iam concentrações variando, aproximada- 
mente, de 625000 a 2000 000 B. coli/ml. 

Nos trabalhos preparatórios para o dimensio- 
namento do grande exutor de Los Angeles, de 
8 km de extensão, medições dos teores de con- 
centração bacteriológica nas várias etapas de 
tratamento de esgoto na estação existente (Hy- 
perion Treatment Plant), indicaram: [10] 


— no esgoto bruto .... 366000 B. coli/ml 
(média geométrica do NMP) 
— ap0s a passagem por 
grelhas e decantadores 259000 Idem, idem. 
— após tratamento bioló- 
gico (lamas activadas) 112000 Idem, idem. 


Gunnerson [11] cita os valores de 500 000 
B. coli/fml e de 1000000 B. coli/ml (médias geo- 
métricas dos NMP), tomados nos cálculos do 
citado exutor, valores que, portanto, pressupõem 
a adopção dum certo factor de segurança em 
relação aos valores medidos e atrás referidos. 

Os caudais de esgoto em Los Angeles são 
muito maiores que os de Luanda; atendendo, 
porém, a que a concentração bacteriológica no 
esgoto varia muito pouco com a capitação [9] 
(varia muito menos do que, por exemplo, o b. 
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O.D., a percentagem de matérias em suspensão, 
etc), mas não deixando de se admitir uma certa 
situação por segurança, tomou-se, no estudo 
feito, o valor de 2000 000 B, coli/ml (média geo- 
métrica dos NMP) para teor da concentração bac- 
teriológica do esgoto bruto. Isto implica, para 
que os teores da concentração bacteriológica nas 
zonas de banhos não excedam 30 B. coli /ml, que 
se garanta uma redução, do teor inicial, igual a 
67 000 vezes. 


3.3 — Redução da concentração bacteriológica. Facto- 
res intervenientes 


A representação, num gráfico, da curva de 
redução da concentração bacteriológica em fun- 
ção do tempo, apresenta, num caso mais geral, 
três fases distintas: [12] 


— uma primeira fase na qual a concentração 
poder-se-á manter constante ou mesmo au- 
mentar ; 

—uma fase seguinte onde se verifica uma 
acentuada redução da concentração ; 

—uma fase final na qual se observa uma 
grande resistência na redução da concen- 
tração. 


A expressão que define a redução da concen- 
tração bacteriológica é do tipo: 


N 


Ú 


— 10—k (t-to) 
na qual: 


No — número de bactérias em t==to (concen- 
tração inicial) ; 

Nt — número de sobreviventes após o tempo t 
(em horas), contado a partir de to (em 


horas); 

to — número de horas correspondentes à pri- 
meira fase; 

k — factor constante 


Estudos realizados com utilização de vários 
tipos de meios ambientes, levaram à conclusão 
de que a primeira fase normalmente existe quando 
a temperatura é baixa, ou quando ao meio se 
adiciona, por exemplo, uma solução de lactose, 
ou ainda quando o meio está sujeito a grande 
agitação; de qualquer modo, a redução de con- 


centração bacteriológica nas águas do mar (e no 

próprio mar, não em ensaios de laboratório), é 

tal que, como regra, na expressão atrás citada 

to==0, ou seja, não existe a primeira fase. 
Assim : 


Ni 

log — 

qu Ns 
ta — tu 


onde: 


Ni4 — número de bactérias no tempo ti; 
Ns — número de bactérias no tempo tz, após 
o tempo lt. 


Se se fizer N, == 10 Ns e ty = 0, virá: 


too 


significando tw, o tempo, em horas, necessário 
para se reduzir de 90% a concentração bacte- 
riológica. 

Em ensaios levados a efeito na Baía de Santa 
Mónica (Los Angeles) concluiu-se que os facto- 
res intervenientes ma redução da concentração 
bacteriológica eram de três ordens: 


— diluição, provocada pela mistura das águas 
do mar e do esgoto nelas lançado; 

— sedimentação das matérias decantáveis cons- 
tituintes do esgoto, na zona em torno da 
saída do exutor; 

— mortalidade, devido à temperatura, acção 
bacteriofágica de certos animais e plantas 
inferiores, etc. 

E desse modo, na expressão: 


N: og 
—— sms TO) 
No ' 
é: 
1 1 1 1 
ad 
to toog too, tum 
sendo: 


tooy — tempo, em horas, para reduzir de 90 “/ 
a concentração bacteriológica por efeito 
exclusivo da diluição; 


toy, — idem, por efeito exclusivo da sedimen- 
tação ; 
too — idem, por efeito exclusivo da morta- 


lidade. [11] 


a) Factor diluição 


As clássicas experiências de Rawn e Pal- 
mer [13] sobre a difusão de jactos líquidos na 
água do mar, levaram os seus autores a propor, 
para jactos de tubagens horizontais, a expressão 
empírica : 
05(L+3)2 


dg Q 05! 


na qual: 


So — diluição, no topo do jacto, na intersecção 
com a superfície do mar, numeéricamente 
igual à relação entre o volume da mis- 
tura da água do mar e líquido diluído, e 
o volume do líquido diluído, numa amos- 
tra (diluição inicial) ; 

Q — cauda do jacto, em galões/min; 

L — desenvolvimento do jacto, segundo o eixo, 
de forma parabólica, em pés. 


Fig. 2 — Difusão de jactos de tubagens horizontais 
(Rawn e Palmer) 


Se se representar por X a distância, medida na 
horizontal, entre a saída do exutor e o ponto de 
intersecção do eixo do jacto com a superfície do 
mar; e se por Y se representar a profundidade, 
em relação à superfície do mar, a que se encon- 
tra a saída do exutor; finalmente, se se fizer 


34 
a= VºD, onde V é o velocidade do jacto à 
saída do exutor em pés/s e D o diâmetro do 
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exutor em pol., poder-se-á escrever, ainda de 
acordo com as mesmas experiências; 


| 2 2 A 
L=Y + (0,8)! — a : quando — > — 
[y X 3 

V 
PNONA E x É 
L=X11+0,9(—| lI..... uando — << — 
[ia (O) Jem do 


No caso estudado em Luanda era sempre, 
como se verá: 


Deste modo, pois, sabidos O, V e D, poder- 
-se-á determinar a diluição do esgoto no topo 
do jacto (diluição inicial). 

Estas experiências foram mais tarde reanali- 
sadas em laboratório pelo próprio Rawn, e ainda 
por Bowerman e por Brooks. [14] Daí se con- 
cluiu ser a diluição inicial do esgoto função da 
profundidade da saída do exutor (Y) e do nú- 
mero de Froude do jacto (F) definido por: 


onde: 


g = A fa (aceleração aparente 


y esgoto da gravidade) 


Os três últimos autores citados não apresen- 
tam qualquer expressão de cálculo da diluição 
inicial, mas apenas uns gráficos de eixos loga- 


, : + À 
ritmos, onde relacionam So com D e com F. 


Se a zona onde descarrega o exutor estiver 
sujeita a fortes correntes, os valores de So dados 
pelos gráficos podem ser muito inferiores aos 
que na realidade se verificam; nestes casos o 
valor médio de So deverá ser determinado 
por: |14| 

Ux bxh 


mem 


Q 


So Sae 


onde: 


U — velocidade da corrente (m's); 

b —largura da mancha de esgoto no topo do 
jacto, medida segundo a perpendicular à 
direcção de U (m); 
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h — altura da mancha do esgoto, na origem 
(m) ; 
OQ — caudal de esgoto (m/s). 


Para além das expressões indicadas, outros 
autores têm proposto várias fórmulas para o 
cálculo da diluição inicial. 

Assim Conway, Pomeroy e Cooley e Harris 
apresentam as seguintes: 


Y2 
So = (Conway) 
Q 
Se =Y (Pomeroy) 
Y 
So = — (Cooley e Harris) 
3D 


onde: 


So— diluição inicial; 

Y — profundidade, em relação ao nível do 
mar, a que se encontra a saída do esgoto, 
em pés; 

O — caudal do jacto, em galões/min; 

D — diâmetro do orifício de saída, em pés. 


Ter-se-á em atenção, em cálculos posteriores, 
cada uma das expressões indicadas. 


x * 


A partir do topo do jacto a mancha de esgoto 
formada terá tendência a deslocar-se de acordo 
com a velocidade de transporte de massa líquida 
dominante (no caso de ausência total de trans- 
porte de massa, a mancha progrediria igual- 
mente segundo todos os raios passando pelo topo 
do jacto). Como consequência, em pontos afas- 
tados do topo do jacto a diluição irá aumen- 
tando. 

Rawn e Palmer, [13,1,10] bem como Conway 
[1,10] por um lado e Pomeroy, por outro, apre- 
sentam expressões, de carácter mais ou menos 
empírico, que permitem o cálculo da diluição 
sequente em função da diluição inicial. Todavia, 
apenas se fará referência e, como consequência, 
apenas se aplicará, a fórmula indicada por Brooks 
[1,10,5], dado que ela foi desenvolvida com bases 
teóricas e experimentais que oferecem bastante 
confiança se bem que com a admissão de certas 
condições simplificadoras. 

Com referência à Fig. 3, partindo da lei da 
difusão linear, a qual traduz que a passagem ou 
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o fluxo do esgoto diluido através da área unitá- 
ria, perpendicular a x, no tempo unitário, é igual 
de 


a—:s , onde < significa o coeficiente de 


X 


. - o Cc ud - Ed 
difusão e E a variação de concentração do es- 
x 


goto segundo x, Brooks admitiu em seguida as 


Fig. 3 — Formação da mancha de esgoto 
à superfície (Brooks) 


seguintes hipóteses, verificadas experimental- 
mente no mar dentro da precisão das medidas: 


—&: é função de L (:=-«Lº, com «=-0,01 em 
unidades c.g.s. segundo Pearson); 

— L é função exclusivamente de x; 

— a diluição segundo o eixo vertical ao plano 
(x, y) é desprezável; 

— o caudal de esgoto é constante. 

Em continuação Brooks deduz: 


Crmáx x o =—— Cs erf 


onde: 


Cmáx « o — Concentração máxima à distância x 
"do topo do jacto (origem); 

Co — concentração no topo do jacto (con- 
centração inicial) ; 

erf — abreviatura, na literatura inglesa, de 

standard error function, representa a 

área, entre — y e +y, da curva 

normal de Gauss relativa à variável 

Y, designando-se por y, e por sim- 

plicidade, o valor da raiz que figura 


na expressão; 
p 72 Es 
É = ; 
Ub 


£o = ab: 


b — largura inicial da mancha de esgoto, numa 
direcção perpendicular a U; 

U — velocidade de transporte da massa lIi- 
quida. 


Deste modo, se, por exemplo, for 1500 000 
B. coli/fml a concentração bacteriológica no esgoto 
bruto e So == 100, a concentração inicial, no topo 
do jacto, será Co; == 15000 B. coli/jml; e sabidos 
b e U, a expressão anterior permitirá o cálculo 
da concentração máxima Cmix, o à uma distância 
x da origem, devida exclusivamente à diluição 
do esgoto. 

De notar que num plano perpendicular ao 
eixo dos xx a concentração apresenta valores 
cada vez mais pequenos nas proximidades dos 
limites da mancha do esgoto. 

Teóricamente, a distribuição da concentração: 
nesse plano genérico, segue a lei de Gauss. 

Consequentemente, a relação entre o valor 
genérico Cx,y desse plano e o valor Cw,o, no eixo 
dos xx, é da forma: [16] 


cado És 
== e 2 o) 


onde 7 representa o desvio-padrão da lei de 
Gauss. 
O seu valor é: 


E 
Exa 


, 


Tabelas existentes [16] permitem a determina- 


a partir de y 
x,0 g 
Como se vê, portanto, a expressão de Brooks 


conduz ao cálculo directo do valor da concentra- 
ção ao longo do eixo dos xx; sabidos estes vá- 
rios valores, poder-se-ão determinar os valores 
da concentração em qualquer ponto x,y da man- 
cha de esgoto, pela aplicação da lei de Gauss. 

Evidentemente, será sempre, nesta aplicação, 


Cx,y 


ção da relação 


1 
y< > em cada plano perpendicular ao eixo 


dos xx. 

O que fica exposto até aqui refere-se a deter- 
minação dos valores da concentração em cada 
ponto da mancha do esgoto; porém, outro as- 
pecto da análise de Brook que interessa focar, 
diz respeito ao cálculo da largura da mancha de 
esgoto ao longo do eixo dos xx, isto é, da 
função L =L (x). 
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Brooks deduz para essa função a seguinte 


forma: 
L ( à à E 
=>=lgí = RE 
b : e 3 b 


na qual já são conhecidos os significados dos 
vários parâmetros, dispensando-se neste lugar a 
sua demonstração. 


b) Factor sedimentação. Factor mortalidade 


Gunnerson, na obra já citada, relata os resul- 
tados de interessantes investigações que realizou 
nas zonas do mar envolvendo as saídas de exu- 
tores que debitavam quer esgoto simplesmente 
decantado como esgoto decantado ou sujeito a 
cloração e ainda esgoto previamente tratado por 
lamas activadas. 

Para esgoto simplesmente decantado, e lançado 
por diferentes exutores em diferentes zonas, de- 
terminou os valores respectivos de ty (devidos 
simultâneamente à diluição, sedimentação e mor- 
talidade), os quais se escalonavam entre 1,4 e 
4 horas. Designadamente, na comparação de dois 
casos distintos, um estudo mais pormenorizado 
indicou: 


to0m too too too 
— esgoto decantado 17,8h 20,0h 53h 34h 


— esgoto de lamas 
activadas ...... 17,8h 20,0 h 21,0h 6,5h 


Estes resultados pôem em evidência o seguinte: 
enquanto os tempos necessários para reduzir de 
90º/) a concentração bacteriológica por efeito ou 
da mortalidade ou da diluição são os mesmos 
nos dois casos, o tempo correspondente ao efeito 
de sedimentação é muito mais pequeno no esgoto 
simplesmente decantado. Aliás é esta uma con- 
clusão absolutamente lógica. 

Se nos valores atrás transcritos se combinarem 
os efeitos de tw, e tm, ter-se-á: 


| 5 1 


ú tuOm too, 


donde, por substituição: 


toom+s * 4 horas no esgoto decantado. 
toom+s = 10 horas no esgoto de lamas activadas. 
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3,4 — Redução da concentração bacteriológica. Fórmu- 
las de cálculo 


Tanto Conway como Pomeroy apresentam fór- 
mulas de cálculo da redução da concentração 
bacteriológica por efeito combinado da diluição, 
sedimentação e mortalidade. Não foram seguidas 
no trabalho a que se refere este artigo dado que 
se considerou a de Brooks de formulação mais 
adequada à realidade. Quanto a Rawn e Palmer, 
na sua já indicada obra, limitam-se a afirmar que 
os teores de concentração bacteriológica admissi- 
veis para zonas de banhos serão atingidos logo 
que a diluição do esgoto iguale 225. Esta indica- 
ção é manifestamente insuficiente e a prática não 
a tem demonstrado. 

À expressão de Brooks referida atrás permite, 
como se disse, o cálculo da redução da concen- 
tração bacteriológica por efeito exclusivo da dilui- 
ção. Todavia, Brooks formulou uma outra hipó- 
tese, a acrescentar às citadas, e que foi a de supor 
que a redução da concentração bacteriológica por 
efeito combinado da mortalidade e sedimentação 
é proporcional à concentração (—k'c). Daí a 
expressão mais geral de Brook, que se escreve: 


Cráxy o = Co e"! erf areas i 
(1 + E Re) — 
N 3 -b 


na qual t=x/U e k = k/log e, dado que: 


O valor de k a tomar para a determinação 
de k' é, como já ficou dito: 


sendo, todavia, neste caso, ty» o tempo, em horas, 
necessário para reduzir de 90º a concentração 
bacteriológica por efeito combinado da mortali- 
dade e da sedimentação. 

Os valores de ty» que se tomaram na aplicação 
da expressão de Brooks ao exutor de Luanda 
estudado, foram fixados por comparação com va- 
lores admitidos nos cálculos de outros exutores 
e com valores medidos nas zonas que circundam 
exutores já construídos, referido na literatura con- 
sultada. Como se viu, o facto do esgoto a lançar 


